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1. Introduction 



One of the most important open problems in hadronic physics (studied since long be- 
fore the discovery of QCD) is to explain/predict the (asymptotic) high-energy behaviour 
of hadron-hadron total cross sections. Present-day experimental observations (up to 
center-of-mass total energy i/i =1.8 TeV) seem to be well described by a pomeron-like 
high-energy behaviour (see, for example, Ref. [1] and references therein): 



This behaviour is known in the literature as the soft pomeron, to be distinguished from 
the well-known BFKL (or hard) pomeron [2], i.e., ~ s~^^°^'^, with = g"^ /An, obtained 
in perturbative QCD. As we believe QCD to be the fundamental theory of strong inter- 
actions, we also expect that it correctly predicts from first principles the behaviour of 
hadronic total cross sections with energy. Anyway, in spite of all the efforts, a satisfactory 
solution to this problem is still lacking. 

We should also remind the reader at this point that the pomeron-like behaviour (II. ip is, 
strictly speaking, theoretically forbidden (at least if considered as a true asymptotic be- 
haviour) by the well-known Froissart-Lukaszuk-Martin (FLM) theorem [3] (see also |1]), 
according to which, for s — * oo, crtot(s) < log^ ( f" ) ' where is the pion mass and Sq 



is an unspecified squared mass scale. In this respect, the pomeron-like behaviour (11. ip can 
at most be regarded as a sort of pre- asymptotic (but not really asymptotic^?) behaviour of 
the high-energy total cross sections (see, e.g., Refs. OEIE] and references therein), valid 
in a certain high-energy range. 

From a general theoretical point of view, the optical theorem (which is a consequence 
of unitarity) allows one to derive hadron-hadron total cross sections from the expressions 
of the hadron-hadron elastic scattering amplitudes M.{hh){s,t) {t being the transferred 
momentum squared): 



High-energy hadron-hadron elastic scattering processes at low transferred momentum 
(usually called soft high-energy scattering processes) possess two different and widely 
separated energy scales: the center-of-mass total energy squared s, which is a hard scale 
(s ^ 1 GeV^: formally we consider the limit s — > oo), and the transferred momentum 
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squared t, which is a soft scale, smaller than (or nearly equal to) the typical energy scale 
of strong interactions (|t| < 1 GeV^ <C s). In this energy regime we cannot fully rely on 
perturbation theory and a genuine nonperturbative treatment is in order. 

Since Nachtmann's seminal paper in 1991 [8] a lot of work has been done on the prob- 
lem of soft high-energy scattering in the framework of nonperturbative QCD P, [101 [HI 
[T2I [13]. Using a functional-integral approach, hadron-hadron elastic scattering ampli- 
tudes are reconstructed from the correlation functions of certain Wilson loops, defined in 
Minkowski space-time (see Section 2 below). As it has been shown in [T^ \TE\ \W \ [T7] (see 
Section 2 below), such correlation functions can be recovered after proper analytic con- 
tinuation from their Euclidean counterparts, i.e., correlation functions of certain Wilson 
loops forming an angle 6 in Euclidean space: this has paved the way to the application of 
nonperturbative techniques, which are normally available only in the Euclidean formula- 
tion of QCD. Some calculations of the loop-loop Euclidean correlation functions already 
exist in the literature (see Section 4 below), using the so-called Stochastic Vacuum Model 
(SVM) [18], the Instanton Liquid Model (ILM) [19], or the AdS/CFT correspondence for 
strongly coupled gauge theories in the limit of a large number of colours [20], [21], [22] . 
The Euclidean correlator is then continued to the corresponding Minkowskian correlation 
function using the above-mentioned analytic continuation in the angular variables. The 
loop-loop correlation functions (both in the Minkowskian and in the Euclidean theories) 
have also been computed exactly in the first two orders of perturbation theory, 0{g'^) and 
C(/), in Ref. [5]. 

Of course, each of these models, used in the calculation of the Euclidean correla- 
tion functions, has its own limitations, which are reflected in the variety of answers in 
the literature: someone finds constant cross sections, someone else finds a soft-pomeron 
behaviour, someone else finds a hard-pomeron behaviour . . . (And maybe the true asymp- 
totic behaviour is log^(s/so), thus saturating the FLM bound!?). Unfortunately, these 
limitations are often out of control, in the sense that one does not know how much infor- 
mation is lost due to the involved approximations. This is surely a crucial point which, in 
our opinion, should be further investigated. (For example, in Ref. [23] it has been shown 
that the Euclidean-to-Minkowskian analytic-continuation approach can, with the inclu- 
sion of some extra, more or less plausible, assumptions, easily reproduce a pre-asymptotic 
pomeron-like behaviour, like the one in (11.11) . which violates the FLM bound.) 

A real breakthrough in this direction could be provided by a direct numerical calcu- 
lation of the loop-loop Euclidean correlation function in Lattice Gauge Theory (LGT). 
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In particular, one is interested in the dependence of the correlation function on the angle 
9 between the loops, from which the energy dependence of the physical scattering am- 
plitudes can be derived after a proper analytic continuation to Minkowski space-time. 
Clearly a lattice approach can at most give the above-mentioned function only for a finite 
set of 6'-values, from which it is clearly impossible (without some extra assumption on 
the interpolating continuous function) to get, by analytic continuation, the correspond- 
ing Minkowskian correlation function (and, from this, the elastic scattering amplitudes 
and the total cross sections). However, the lattice approach, which provides a real first- 
principle determination of the loop-loop correlator in Euclidean space, can be used to in- 
vestigate the goodness of a given existing analytic model (such as SVM, ILM, AdS/CFT, 
and so on . . .) or even to open the way to some new model, simply by trying to fit the 
lattice data with the considered model. This is exactly what we shall try to do in this 
paper. 

The plan of this paper is the following. In Section 2 we give (for the benefit of the 
reader) a brief review of how high-energy scattering amplitudes can be reconstructed, us- 
ing a functional-integral approach, in terms of certain correlation functions of two Wilson 
loops and we also briefly recall some relevant analyticity and crossing-symmetry proper- 
ties of these loop-loop correlation functions, when going from Euclidean to Minkowskian 
theory. In Section 3 we shall see how these (Euclidean) loop-loop correlation functions 
can be evaluated in lattice QCD and discuss the main technical complications. In Section 
4 we shall compare our numerical results with some analytical nonperturbative estimates 
that appeared in the hterature, discussing in particular the question of the analytic con- 
tinuation from Euclidean to Minkowskian theory and its relation to the still unsolved 
problem of the asymptotic s-dependence of the hadron-hadron total cross sections. In 
Section 5 we draw our conclusions and show some prospects for the future. 



2. Hadron-hadron (dipole— dipole) scattering 
amplitudes from Wilson loop correlators 

It was shown in Refs. P, [THl [H] (for a review see Refs. [H [12]) that the high-energy 
meson-meson elastic scattering amplitude can be approximately reconstructed in two 
steps. 
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i) One first evaluates, in the functional-integral approach, the high-energy elastic 
scattering amplitude M.(dd) of two qq pairs (usually called dipoles) of given transverse 
sizes Ri± and R2± and given longitudinal-momentum fractions /i and /2 of the two 
quarks in the two dipoles respectively. This dipole-dipole scattering amplitude turns out 
to be governed by the (properly normalised) correlation function of two Wilson loops Wi 
and VV2, which follow the classical straight lines for quark (antiquark) trajectories: 



M(dd){s,t; Ri_lJi, R2±, f2) = -i2s d z±e 



{m){m) 



(2.1) 



where s and t = —\qj_\'^ {qi_ being the transferred momentum) are the usual Mandelstam 
variables. The expectation values (W1W2), (Wi), (W2) are averages in the sense of the 
QCD functional integral {{0[A]) = ^ J[dA]0[A]detQ[A]e'^^^\ Z = J[dA]detQ[A]e''^^'^\ 
with S[A] the pure-gauge action and detQ[74] the fermion-matrix determinant) and the 
Wilson loops Wi and W2 are defined as follows (for Ai'c-colour QCD): 



— Tr < Pexp 



—ig f A^{x)dx'' 



Cl,2 



[2.2) 



where V denotes the path ordering along the given path C and A^ = A^T", T"" being 
the generators of the SU{Nc) Lie algebra in the fundamental representation; Ci and C2 
are two rectangular paths which follow the classical straight lines for the quark [Xg(r), 
forward in proper time r] and the antiquark [Xq(r), backward in r] trajectories, i.e., 



Ci 
C2 



z 



m 



X. 



2q 



T 



m 



r + (l-/2 



X2',{r) 



m 



T 



f2R^, 



(2.3) 



and are closed by straight-line paths at proper times r = ±T, where T plays the role of 
an infrared (IR) cutoff [16], which can and must be removed in the end, by letting 
T ^ 00. (In fact, differently from the parton-parton scattering amplitudes, which are 
known to be affected by IR divergences, the elastic scattering amplitude of two colourless 
states in gauge theories, e.g., two qq meson states, is expected to be an IR-finite physical 
quantity j25j.) 

Here pi and p2 are the four-momenta of the two dipoles, taken for simplicity with the 
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same mass m, moving (in the center-of-mass frame) with speed V and —V along, for 
example, the x^-direction: 

= 772 (^cosh^,sinh|,Oj_j , pa = (cosh |, - sinh |, 0_l) • (2.4) 

Here X = 2 arctanhl^ is the hyperbolic angle between the two trajectories Iq and 2q, i.e.. 
Pi ■ P2 = fTT-"^ coshx. In the high-energy limit s ^ oo (i.e., x ~^ +oo)) 

s = (pi +p2)^ = 2m^ (coshx + 1) , i.e.: x ~ log ( — ^ ) • (2-5) 

Moreover, Ri = (0,0, ^i_l), R2 = (0,0, ^2±) and z = (0,0, ^_l), where z± = {z'^,z^) is 
the impact-parameter distance between the two loops in the transverse plane. The two 
Wilson loops are schematically shown in Fig. [H 

If we introduce the following notation for the normalised connectec? loop-loop correlator 
(in the presence of a finite IR cutoff T) , Q 

where the arguments "1" and "2" stand for "-Ri_l,/i" and "-R2±,/2" respectively, the 
dipole-dipole scattering amphtude (12. ip can be rewritten as 

M(dd)is,t; Ri^, fi, R2±, f2) = 

-I 2s ! d^z^e^'^^-'^Guix ~ log('^);T^+oo;il;l,2). (2.7) 

ii) The hadron-hadron (in our case meson-meson) elastic scattering amplitude A4(hh) 
can then be obtained by averaging the dipole-dipole scattering amplitude (12.71) over all 
possible dipole transverse separations Ri± and R2± and longitudinal-momentum fractions 
/i and f2 with two proper squared hadron wave functions \'ipi{Ri±, fi)\^ and \'ip2{R2±, /2)P, 
describing the two interacting mesons: 

M {hh){s,t) = J d^Ri± J dfi \ipi{Ri±, fi)\'^ J d'^R2± J df2 \4'2{R2±, 12)1"^ 

X M^dd){s,t;Ri^,f,,R2^,f2). (2.8) 



"The quantity Gm in Eq. (US]) is equal to the quantity Gm - 1 of Refs. pJl fTTl [23 l [26] . 
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(For the treatment of baryons, a similar, but, of course, more involved, picture can be 
adopted, using a genuine three-body configuration or, alternatively and even more simply, 
a quark-diquark configuration: we refer the interested reader to the above-mentioned 
original references P El [IDl [IH [12] •) 

All the above refers to the theory in the (physical) Minkowski space-time. However, 
as we have already observed in the Introduction, most of the nonperturbative methods in 
field theory are available only in the functional-integral approach formulated in Euclidean 
space: therefore we are interested in the problem of the analytic continuation of the loop- 
loop correlation functions into/from Euclidean space. In the Euclidean theory we can 
consider the correlation function of two Euclidean Wilson loops Wi and W2, defined in 
the same way as they are in Minkowski space-time, taking into account that the metric 
is now Euclidean rather than Minkowskian, and running along two rectangular paths Ci 
and C2 which follow the following straight-line trajectories 

C, : X|;^^(r) = ^r+(l-/2)i?2i.„ X^^(r) = - /^ii^i.,, (2.9) 

and are closed by straight-line paths at proper times r = ±T. Here the Euclidean 
coordinates are Xe^ = {Xei^ Xe2, Xe^,, Xea), where Xea corresponds to the "Euclidean 
time", and -Ri^; = (0,i?i_L,0), R2E = (0,-R2±,0) and ze = (0, ij_,0), where -Ri_l, R2± and 
z± are exactly the same transverse vectors introduced in the Minkowskian case. Moreover, 
in the Euclidean theory we choose the four-vectors piE and p2E to be 

/ e ^ e\ f e ^ e\ , , 

Pie = "2 I sm-,0^,cos- I , p2i? = I - sm -, 0_l, cos - 1 , (2.10) 

6 being the angle formed by the two trajectories Iq and 2q in Euclidean four-space, i.e.. 
Pie ■ P2E = cos 0. 

It turns out [161 [IT] that the Minkowskian quantity Qm [with x ^ l^"*"] can be reconstructed 
from the corresponding Euclidean quantity [with 9 G (0, vr)] 

(where (. . .) are now averages in the sense of the Euclidean functional integral, i.e., 
{0[Ae\) = ^ J[dAE]0[AE]detQE[AE]e-'^^^^\ Ze = J[dAE]detQE[AE]e-''^^^^\ with 
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S'£;[A£;] the pure-gauge Euclidean action and detQcfAe] the Euchdean fermion-matrix 
determinant jll by an analytic continuation in the angular variables 6 —ix and in the 
IR cutoff T — > iT, exactly as in the case of Wilson lines [TH [151 [IS]. This result [161 [H] 
is derived under certain hypotheses of analyticity in the angular variables and in the 
IR cutoff T. In particular, one makes the assumption [26] that the function Qe, as a 
function of the complex variable 6, can be analytically extended from the real segment 
(0 < Re^ < TT, Im^ = 0) to a domain Ve, which also includes the negative imaginary 
axis (Re6' = 0+, Im6' < 0); therefore, the function Qm, as a function of the complex vari- 
able can be analytically extended from the positive real axis (Rex > 0, Imx = 0+) 
to a domain Dm = {x ^ \ —^X^^-e}, which also includes the imaginary segment 
(Rex = 0, < Imx < vr). The validity of this assumption is confirmed by explicit cal- 
culations in perturbation theory [T^ [T71 [5]. The domains Ve and Vm are schematically 
shown in Fig. [2l Denoting with Qm and Qe such analytic extensions, we then have the 
following analytic-continuation relations [T71 [26] : 

^e(^;T;?^;1,2) = ^m(^^; -^T; ^1; 1, 2), G P^,; 

gMix;T;zrA,2) = ^s(-zx;^T;il;l,2), Vx G Pm- (2.12) 

As we have said above, the loop-loop correlation functions (both in the Minkowskian and 
in the Euclidean theories) are expected to be IK-finite quantities, i.e., to have finite limits 
when T —>■ oo, differently from what happens in the case of Wilson lines. One can then 
define the following loop-loop correlation functions with the IR cutoff removed: 

Cm(x;^±;1,2) = lim ^M(x;7';il; 1,2), 

CEi9;zr,l,2) = \im gE{9;T;z^;l,2). (2.13) 

It has been proved in Ref. [IT] that, under certain analyticity conditions in the complex 
variable T [conditions which are also sufficient to make the relations fl2.12p meaningful], 
the two quantities fl2.13p . obtained after the removal of the IR cutoff (T oo), are still 

^We note here that, exploiting the invariancc of the theory under charge conjugation (C), one can show 
that the correlation function Ge, 'd^s well as the quantity Ce to be defined below in Eq. (|2.13p . is a real 
quantity (differently from the corresponding Minkowskian correlator). In fact, under a C transformation, 
AEfi j4^^ = —A*-^^ ~ ^^Efi^ the integration measure [dA^;], the gauge-field action Se and the fermion 
matrix determinant detQ£;[A£;] (which arc real quantities) are invariant, while a Wilson loop W goes into 
its complex conjugate W*. 
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connected by the usual analytic continuation in the angular variables only: 



Cm(x;^±;1,2) = CE{-tx;zrA,2), Vx g Pm- (2.14) 

This is a highly non-trivial result, whose general validity is discussed in Ref. [17]. The 
validity of the relation (12.141) for the loop-loop correlators in QCD has been also verified in 
Ref. [S] by an explicit calculation up to the order 0{g^) in perturbation theory. However 
we want to stress that the analytic continuation (12.121) or (I2.14p is expected to be an 
exact result, i.e., not restricted to some order in perturbation theory or to some other 
approximation, and is valid both for the Abelian and the non-Abelian cases. 

It has also been recently shown in Refs. [261123] that the analytic-continuation relations 
( ]2.12p allow us to deduce non-trivial properties of the Euclidean correlator Qe under the 
exchange 9 ^ ti — 9 and of the Minkowskian correlator Qm under the exchange x ~^ i'^~Xi 
corresponding to the exchange from a loop-loop correlator to a loop-antiloop correlator, 
where an antiloop is obtained from a given loop by exchanging the quark and the antiquark 
trajectories: 

GEiTT - 9; T; zr, 1,2) = gE{9; T; z^, 1,2) = 6;^(^; T; zl; T, 2), \/9 E (0, vr); 

GMiiTi - x; T; zr, 1, 2) = QMix; T- z^- 1, 2) = 6?m(x; T\ T, 2), Vx e M+, (2.15) 

where the arguments "1" and "2" stand for "— -Ri±, 1 — /i" and "— -R2±, 1 — /2" respectively. 
These two relations are known as crossing-symmetry relations for loop-loop correlators. 
As they are valid for every value of the IR cutoff T, completely analogous relations 
also hold for the loop-loop correlation functions Cm and Ce with the IR cutoff removed 
(T ^ cx)), defined in Eq. fl27[3|) : 

Ce{7i - 9; zr, I, 2) =CE{9-zr, 1,2) =CE{9;zr, 1,2), V^G (0,7r); 

CM{t7i - x; ZL, 1, 2) = Cm[x; zr, 1,2) = Cm{x; z±; T, 2), Vx g m+. (2.I6) 

Taking into account the aforementioned analytic-continuation relations, we can rewrite 
Eq. (12. 8p in terms of the Euclidean correlation function C^; as 

Mi^Hh){s,t) = -12S j d'Rl±j^ dh \MRl±Jl)\' j ^'^2±^ df2 \MR2±,f2)\' 

X J d^z^e''^^-'^CE{9^-i\og(^^y,zr,Ri±,fi,R2±,f2). (2.17) 

9 



By virtue of the optical theorem, Eq. fll.2|] . the total cross section is then given by the 
expression 



aS\s) ~ -2 / d'R,^ [ dh \MRi±Ji)\' I d'R2± [ dh \MR2LJ2)? 
J Jo J Jo 

X /rfVReC,,(^^-zlog(^);^x;i?i±,/i,i?2±,/2). (2.18) 

If one chooses hadron wave functions invariant under rotations and under the exchange 
fi—*l — fi (see Refs. [121 [13] cind also [1], §8.6, and references therein), the correlation 
function Ce in Eqs. f l2.17p and fl2.18p can be substituted (without changing the result) 
with the following averaged correlation function: 

C^*'^(6'; /i, |^2±|, /2) = J dRi_L J dR2± 

X- |C£;(6'; z±; Ri±, fi, R2±, /2) + Ce{0; ^1; Ri±, 1 — /i, R2±, /2) 

+ CEie; il; ^1^, A, ^2±, 1 - /2) + CEie; z^- ^1^, 1 - A, ^2±, 1 - A) } , (2.19) 

where / dRi± stands for integration over the orientations of Ri±. We note here that, as a 
consequence of the crossing-symmetry relations Eq. (]2.16p . the function is automat- 
ically crossing-symmetric, i.e., C^^iji — 9; . . .) = C'^^{6] . . .) for fixed values of the other 
variables. 

In the following, if not specified otherwise, we will take for simplicity the longitudinal- 
momentum fractions /i and A of the two quarks in the two dipoles (and, therefore, also 
the longitudinal-momentum fractions 1 — /i and 1 — A of the two antiquarks in the two 
dipoles) to be fixed to 1/2: this is known to be a good approximation for hadron-hadron 
interactions (see Refs. [H [T2] and references therein). We will also adopt the notation 
Qe{0] T; zr, Ri±, R2±) = Ge{0] T; z_i; Ri^, fi = |, R2±, A = |), and similarly for Ce and 
C'^'^. Note that in this case the quantity C^^ reduces to the average over the transverse 
orientations only: 

Cr(^;^±;|^i±|,|^2±|) = J dRi^ J dR2±CEie;zr,Ri±,R2i.). (2.20) 

In the next section we will show how it is possible to calculate these Euclidean correlation 
functions in LGT and discuss the main technical difficulties of this approach. 
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3. Loop— loop correlators on the lattice 



Being a gauge-invariant quantity, the Wilson-loop correlation function is a natural 
candidate for a lattice computation; nevertheless, a number of complications arise because 
of the explicit breaking of 0(4) invariance on a lattice. The major complication is due to 
the limited number of possible orientations of rectangular loops on the lattice. As straight 
lines on a hypercubic lattice can be either parallel or orthogonal, the values of 9 directly 
accessible are limited to 0°, 90° and 180°. To cover a significantly large set of angles, we 
then have to make use of off-axis (and so non planar) Wilson loops, thus introducing in 
our approach another approximation that should be carefully discussed. 

The loops involved in the calculation of have one side (from now on the longitudinal 
side) in the {xei.xea) (longitudinal) plane and the other {transverse side) in the [xe2^ xez) 
(transverse) plane and their centers are separated in the transverse plane, so that the 
problem of reproducing the loop configuration consists effectively of two distinct two- 
dimensional problems. It is reasonable to evaluate the loop sides on the lattice paths 
that minimise the distance from the true, continuum paths, in order to stay as "close" 
as possible to the continuum limit, adopting essentially the same strategy as in computer 
graphics when drawing straight lines on a screen (see Fig. [3]); such ^^minimal-distance 
paths" can be found in a very efficient way by means of the well-known Bresenham 
algorithm [27], which has already been used in lattice calculations (see e.g. [28], where 
it is also generalised to the three-dimensional case). To every continuum straight path 
having as endpoints two coplanar lattice points we can then associate unambiguously a 
Wilson line on the lattice by means of the minimal-distance prescriptiorQ and then build 
with them the Wilson loops we are interested in. Such loops are identified by the position 
of the center and by two two-dimensional lattice vectors; we thus define Wl(/||; rj_; n) to be 
the lattice Wilson loop evaluated on the minimal path that approximates the rectangle 
having as corners the lattice points n — 1/2 — r/2, n + 1/2 — r/2, n + 1/2 + r/2 and 
n — 1/2 + r/2, where n, / and r are vectors in lattice units, n is the position of the center 
ancill /|| = (/i, ^4), 12 = 13 = 0, fj_ = (r2, rg), n = r^ = 0. 

*For lines inclined at 45° with respect to an axis a certain ambiguity remains, but we can average over 
equivalent paths. 

''^The components of I and r must be integers. Note that the center of the loop may not lie on a lattice 
point. 
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On the lattice we then define the correlator 



rn r X - - ^ (W^(ri||;fi^;c;)W^(r2||;r2x;0)) . 

where d = {0,dj_,0), d± = (^2,(^3). As (full) 0(4) invariance is broken, this correlator 
depends explicitly on the (two-dimensional) lattice vectors fij_ {i = 1, 2) and dj_ rather 
than on their scalar products: indeed, for a given relative orientation we can find different 
realisations in terms of lattice vectors, generally inequivalent as they involve different 
Wilson-loop operators. Anyway, as rotation invariance is restored in the continuum limit, 
we expect 

GL{h\\j2\\;d±;ri±,f2_i) ^ gE{0;Ti = aLi/2,T2 = aL2/2;ad±;afi±,af2±), (3.2) 

where Li = are what we define to be the lengths of the longitudinal sides of the loops 
in lattice units (from now on, ^^lengths"), and ■ l2\\ = L1L2 cos 6*. As it is impossible, in 
general, to have Li = L2, we have relaxed this condition and have considered in Eq. (13. 2p 
the correlation function with two IR cutoffs Ti and T2, 

QEiO; Ti, T2; zr, Ri±, R2±) = ^—(To^^^p)^ - 1' (3.3) 

where with a small abuse of notation we have kept the same notation for the correlation 
function as in Eq. (12. lip . As already pointed out in Section 2, the correlator Qe^^ expected 
to be an IR-finite quantity, so that we can regularise it considering loops with different 
lengths Ti and T2 and taking the limits Ti,T2 00 independently, obtaining the same 
function Ce defined in Eq. (12.131) . 

On top of these field theoretical complications we have to face the numerical difficulty 
(or even feasibility) of a lattice "measurement" of the relevant correlation function of 
Wilson loops. The Wilson-loop expectation value is known to obey an "area law", i.e., 
to vanish exponentially with its area for large areas; moreover, in the 't Hooft large- iVc 
expansion with g'^Nc kept constant such correlators are known to factorise to leading 
order, thus giving 

, _ (vvr'ivf ') _ _ {W}{wP} + o(^) _ j_ 
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We then expect to be a small quantity, obtained from the ratio of two exponentially 
decreasing quantities as T becomes large. Moreover, general considerations make plausible 
a behaviour of the kind Qe ^ e"''''^^' at large separations, and thus we expect the noise 
to overcome the signal after a few lattice spacings. 

We can try to maximise the information obtained from each thermalised configuration, 
in order to reduce the statistical noise, by exploiting the symmetries of the lattice. It is 
easy to see that the chosen prescription for the construction of loops is consistent with 
lattice rotations and reflections [i.e., the cubic subgroup of 0(4)], in the following sense. 
One can perform a cubic transformation on a Wilson loop (in the continuum), and then 
construct its lattice approximation with the given prescription, or alternatively construct 
first the lattice approximation and then perform the same cubic transformation (this 
time on the lattice), and in both cases one would obtain the same result. We can then 
average over cubic transformations of the whole loop configuration and, imposing periodic 
boundary conditions, we can average over lattice translations as well; to further clarify 
the numerical signal we can also use the identities 

GL{h\\,k\\]d±;ri±,r2±) = GL{h\\, k\\] -dr, - fij_, - f2±) 

= GL{h\\,k\\] -dj_;rij_,r2i.) = GL{h\\,k\\;dr, -rij_, -r2±), (3.5) 

which can be proved using invariance under cubic symmetry and the trivial fact that 

Wl(/1; -r±; n) = WLH\\;r±; n). (3.6) 



4. Numerical results 



We have performed a Monte Carlo calculation of the correlation function Ql oi two 
Wilson loops for several values of the relative angle, various lengths and different con- 
figurations in the transverse plane, on a 16^ hypercubic lattice with periodic boundary 
conditions. The link configurations were generated by means of a mixture of (pseudo)- 
heatbath [29l [30l [31] and overrelaxation steps [32] with the usual Wilson action for SU (3) 
pure-gauge theory [33], also known in the literature as the quenched approximation of 
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QCD, which consists in neglecting dynamical fermion loops by setting the fermion matrix 
determinant to a constant. Even though, of course, we cannot exclude that the inclusion 
of dynamical-fermion effects, via the fermion-matrix determinant, could introduce new 
features in the data, in this paper (which, we want to stress, is the first to approach this 
problem from the point of view of lattice QCD) we have preferred to make the easiest 
and most convenient choice, that gives us the possibility of collecting large statistics (also 
considering the various difficulties in measuring the correlation functions Ql, as explained 
in the previous section). 

We have measured the correlation functions (WliWl2) and the loop expectation val- 
ues (VVli) and (Wls), with Wli = yi^L{h\\;ri±;d) and Wl2 = y^Lik\\;r2±;0), on 30000 
thermalised configurations at (3 = = 6.0. As it is well known, the lattice spacing a 
is related to the bare coupling constant g (i.e., to j3) through the renormalisation group 
equation. The lattice scale, i.e., the value of a in physical units, is determined from the 
physical value of some relevant (dimensionful) observable like the string tension or the 
static qq force at some fixed distance {^^Sommer scale") (see e.g. Ref. [M] and references 
therein): in our case one finds that a{(3 = 6.0) ~ 0.1 fm. The choice of /3 = 6.0 on a 16'* 
lattice is made in order to stay within the so-called scaling window": in this sense we 
are relying in an indirect way on the validity of the relation (13.21) between Wilson-loop 
correlation functions on the lattice and in the continuum (and therefore we shall use the 
notation Qe/Ce of the continuum in all the figures reporting our lattice data). An explicit 
test of scaling in our case is more difficult, as one has to keep a large number of length 
scales under control while varying the lattice spacing. A possibility could be to integrate 
over the distance and the sizes of the loops and then study the scaling properties of the 
resulting quantity, although this seems to be a very hard task. 

To keep the corrections due to 0(4) invariance breaking as small as possible, we have 
kept one of the two loops on-axis and we have only tilted the other one as shown in Fig. IH 
the on-axis loop Wli is taken to be parallel to the xei axis, = (Li,0), and of length 
Li = 6,8. We have used two sets of off-axis loops Wl2 tilted at tan"^(l/2) ~ 26.565° 

— * 

and 45° with respect to one of the longitudinal axes; the corresponding lattice vectors l2\\ 
are listed in Table [T], together with their length and the angle formed with /ly. We have 
used loops with transverse size \rij_\ = |r2±| = 1 in lattice units; the loop configurations 
in the transverse plane are those illustrated in Fig. [5l namely d± \\ rij_ || f2± (which we 
call "zz2") and d± ± ri± \\ r2± ("^?/?/")- explained in Section 2, it is interesting to 
also measure the orientation-averaged quantity ("ai^e") defined in Eq. (12.201) . The lattice 
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version of this equation is easily recovered for even (integer) values of the transverse sizes; 
in our particular case, |rij_| = 1, we have to use a sort of smearincj' procedure, averaging 
nearby loops as depicted in Fig. [51 Note that in doing so we are actually averaging over 
the orientations and over the values /j = 0, 1 (z = 1, 2) of the longitudinal-momentum 
fractions, according to Eq. (12.191) . 

As explained in Section 2, we are interested in the T — oo limit and so we have 
to somehow perform it on the lattice. In practice, we have to look for a plateau of the 
correlation function plotted against the loop lengths Li and L2: in Fig. E] we show the 
dependence of the correlator on the length Li = L2 = L of the loops at ^ = 90°. Of 
course, on a 16^ lattice it is difficult to have a sufficiently long loop while at the same 
time avoiding finite size effects and at best we can push the calculation up to L = 8; 
nevertheless, a plateau seems to have been practically reached at about L = Lpi ~ 6-8. 
As 9 varies from 90° towards 0° or 180°, we expect Lpi to gro"ft0; however, the plots in 
Figs. U\M show that the correlation function is already quite stable against variations of 
the loop lengths at Li,L2 — 8 (at least for 6 not too close to 0° or 180°) and so we can 
take the data for the largest loops available as a reasonable approximation of Cl, defined 
as the asymptotic value of as Li, L2 — * 00. We estimate the uncertainty in Cl due to 
this approximation from the variation of Ql with the lenghts Li, L2 as (using the notation 
introduced in Table [1]) 

{SCL)syst = ^{\QL{Li = 8,set = 2) - gL{Li = 8, set = l)\ (4.1) 
+ |e^L(^i = 8,set = 2) - gL{Li = 6,set = 2)\}. 

The errors shown in Figs. [7H9] are the statistical ones only; in all the other figures we 
show the total error obtained adding (in quadrature) the statistical and the "systematic" 
errors defined above. 

As already noticed above, the data are less stable for 9 near 0° and 180°: this is a 
consequence of the relation between the correlation function Qe and the static dipole- 
dipole potential V12 [35] . 



0]T;z±]Ri±,R2±) exp 

T — >oo 



-2TV^2{z±.Ri±.R2 



- 1 (4.2) 



(note that we have set again Ti = T2 = T), from which we expect Qe diverge aX 6 = 0°; 



*Indeed, Lp\ blows up at 0°, 180° due to the relation between the correlation function and the static 
dipole-dipole potential to be discussed below, see Eq. (|4.2p . 
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26.565° 


45° 


63.435° 


90° 


116.565° 


135° 


153.435° 


hw set 1 


(4,2) 


(4,4) 


(2,4) 


(0,6) 


(-2,4) 


(-4,4) 


(-4,2) 


set 2 


(8,4) 


(6,6) 


(4,8) 


(0,8) 


(-4,8) 


(-6,6) 


(-8,4) 


L2 set 1 


2^5 


4v^ 


2v^ 


6 




4V2 


2^5 


set 2 


AVE 


qV2 


aV5 


8 


4V5 


6V2 


4V5 



Table 1: Longitudinal vector /211 and length L2 for the various angles and for the two sets 
of off-axis loops Wl2- 

by virtue of the crossing-symmetry relations f l2.15p . a similar singularity is also expected 
at = 180° [26]. In the following we will consider only 6 ^ 0°, 180°. 

We have considered the values d = 0,1,2 for the distance between the centers of the 
loops. As expected (see the previous Section), the correlation functions vanish rapidly as 
d increases, as can be seen in Fig. [TOl thus making the calculation with our simple "brute 
force" approach very difficult at larger distances. 

From now on we will discuss the issue of the angular dependence of the correlation 
function. As already pointed out in the Introduction, numerical simulations of LGT 
can provide the Euclidean correlation function only for a finite set of 6'-values, and so 
its analytic properties cannot be directly attained; nevertheless, they are first-principles 
calculations that give us (inside the errors) the true QCD expectation for this quantity. 
Approximate analytic calculations of this same function then have to be compared with the 
lattice data, in order to test the goodness of the approximations involved. The Euclidean 
correlation functions we are interested in have been evaluated in the SVM [18], in the 
ILM [19] and using the AdS/CFT correspondence (in the "conformal" [20] and "non- 
conformal" cases [21]). The SVM gives a well-defined quantitative prediction, that can 
be tested numerically against our data; in the ILM and conformal-AdS/CFT cases we 
have qualitative knowledge of the functional dependence of the correlation functions on 
the angle 6 and we can test at least the goodness of the functional form with a fit to the 
lattice data; in the non-conformal-AdS / CFT case the explicit ^-dependence is unknown 
and we are unable to make any comparison with our data. 

In the SVM [18] the Wilson-loop correlation function in iVc-colour QCD [N^ = 3 in 
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our case) is given by the expressioiJll 




A^svM cot e 



KsYM cot 6 



1 



(4.3) 



where -R'svm is a function of z±, Ri± and i?2± only, whose precise expression, which 



we have used to numerically evaluate the correlator (14. 3p in the relevant cases, is given 
in Ref. [18]. The comparison of the SVM prediction for Ce with our data is shown in 
Figs. [TTHT^ in the "ave" case the comparison is made with C"'^'^ as defined in Eq. (12.191) . 
with /i = /2 = (see discussion above). Although in some cases the agreement is quite 
good, at least in the shape of the curve and in the order of magnitude, in general it is far 
from being satisfactory; in particular, in the "21/1/' case for d ^ the SVM prediction is 
orders of magnitude smaller than the lattice results. One can also go the other way round, 
namely try to determine -R'svm in the "2:^2" and ^^zyy" cases with a one-parameter best-fit 
to the data: the results are shown in Figs. [TTHrSl In general, the difference between the 
predicted and the fitted values for i^svM is positive and in some cases the discrepancy 
is larger than 20%, which can be taken as the accuracy level of the model parameters 
(see Ref. [IE])- We note that for c? = 0, 1 the main contribution to the value of -ft'svM, 
evaluated using the SVM expression given in Ref. [18] (which, we recall, consists of a 
perturbative plus a nonperturbative component) comes from perturbative effects, while at 
d = 2 nonperturbative effects are equal to or greater than the perturbative ones. We have 
also tried a best-fit with the following simple functional form: 



which is exactly what one obtains in leading-order perturbation theory [5l [T71 [18] . Notice, 
however, that the coefficient i^pert in Eq. (14.41) can also receive nonperturbative contribu- 
tions, as one can see, for example, when expanding the exponentials in the SVM expression 
n to first ordei§. The values of the chi-squared per degree of freedom {x\ o f ) the 
various fits that we have performed are listed in Table [21 

^Here and in the following formulae we omit the variables other than on which Ce depends. 
■i-This makes sense when i^svM is small, as it happens for example in the large-iVc expansion, where 

it is if SVM = 0(1^ ). 



(4.4) 
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One-instanton effects in the ILM give the following analytic expression for the corre- 
lation function [19]: 

C|^^)(^) = (4.5) 

^ ^ ^ sin^ ^ ^ 

Lattice data are not well fitted by such a function, at least at = 0, 1; at c? = 2 practically 

all the fits we have tried are good, as one can see in Table but we interpret this as 

the result of the large data errors. We can largely improve the fits by adding a term 

proportional to (cot 6')^, 



cS'^''''\0) = ^ + K'{cot9r, (4.6) 



sine ^^^ILMpl 

which describes contributions both from two-instanton [in] and, as we have already said 
above, from leading-order perturbative effects; in the following, we will refer to Eq. (14. 6 p as 
the "ILMp" expression. Differently from the S VM case (and also from the AdS / CFT case 
to be discussed below), the dependence on 9 is not affected by the average over transverse 
orientations (and longitudinal-momentum fractions), so it makes sense to try to also fit 
the "averaged" data with these functional forms. The resulting best-fit functions in the 
d = 0,1 cases are plotted in Figs. [H] and [13 As the fit gives -f^jLMp ^ -^ilmp' seems 
more likely that -R'i'lmp is dominated by perturbative effects; while -R'i'lmp is an order of 
magnitude larger than -f^iLMp at ci = 0, suggesting that perturbative effects are dominant, 
at d = 1 the two parameters are comparable, and at ci = 2 Kilmp becomes larger than 
-^iLMp' suggesting that we are entering the nonperturbative region. 

Using the AdS/CFT correspondence, one obtains for the J\f = 4 SYM theory at large 
Nc and large 't Hooft coupling, and at large distances between the loops [20] 

^^Ads/CFT)^^^ =exp|i^i-i- + K2Cot^ + K3Cosecote| -1. (4.7) 

I sill u J 

Although, of course, there is apparently no reason for this formula to apply in our case, 
we nevertheless have tried to fit the lattice data (in the and ^^zyy" cases) also with 

this functional form; the results are shown in Fig. [161 The fit is not good at (i = 0; it 
may seem good aX d = 1, 2, but we must remember that this is a three-parameter fit, i.e., 
there are only four degrees of freedom: the ILMp best-fit, which has only two parameters 
(i.e., five degrees of freedom) achieves a better or comparable xlof ^^^^ moreover, we 
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d = 







d=l 






d = 2 






zzz/ zyy 


ave 


zzz 


zyy 


ave 


zzz 


zyy 


ave 


SVM 


51 




16 


12 




1.5 


2.2 




pert 


53 


34 


16 


13 


13 


1.5 


2.2 


4.5 


ILM 


114 


94 


14 


15 


45 


0.45 


0.35 


1.45 


ILMp 


20 


9.4 


0.54 


0.92 


1.8 


0.13 


0.12 


0.19 


AdS/CFT 


40 




1 


0.63 




0.14 


0.065 





Table 2: Chi-squared per degree of freedom for a best-fit with the indicated function. 

have explicitly tested some modified versions of the SVM expression (14.31) with two and 
three parameters in the exponents, obtaining smaller values for Xd o f • 

As we have said in the Introduction, the main motivation in studying soft high-energy 
scattering is that it can lead to a resolution of the total cross section puzzle, so it is worth 
discussing what the various models have to say on this point. Using Eq. fl2.18p it is easy 
to see that the SVM, the ILM and the lowest order perturbative expressions give con- 
stant cross sections at high energjii, as in these cases the high-energy limit can be carried 
over under the integral sign, so that the knowledge of the ^-dependence of the corre- 
lation function is sufficient to completely determine the high-energy behaviour of total 
cross sections. However, it may happen that the integrand in Eq. fl2.18p has no definite 
asymptotic behaviour at large s, as it happens, e.g., for the AdS/CFT expression (14.71) . so 
that the remaining integrations have to be carried out before taking s ^ oo. Anyway, as 
the experimentally observed universality suggests that the high-energy behaviour is not 
affected by the details of the hadron wave functions (and thus by the detailed dependence 
of the correlation function on and /j), we would expect that the high-energy asymp- 
totics could be read off after integrating only over the distance \z_[_\ between the loops 
and this would require the detailed knowledge of the dependence of the correlators on 
\z±_\. As already noticed above, however, lattice data show that the relevant correlation 
functions vanish rapidly with the distance: due to our limited numerical knowledge we 
postpone the detailed discussion of the dependence on the distance and related issues to 
subsequent works. 

As a final and important remark, we note that our data show a clear signal of C-odd 
contributions in dipole-dipole scattering. As shown in Ref. [26] (and briefiy recalled in 

§ Actually, the ILM expression ()4.5p for the correlation function (where we have to recall that Kylm is 
a real function) results in an exactly zero total cross section, when inserted in Eq. (|2.18p . 
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Section 2) the loop-antiloop correlator at angle 9 in the Euclidean theory (or at hyperbolic 
angle x iii the Minkowskian theory) can be derived from the corresponding loop-loop 
correlator by the substitution 6 ^ n — 6 (or x — > ztt — x in the Minkowskian theory). 
Because of this crossing-symmetry relation, it is natural (see Ref. [23j) to decompose the 
Euclidean correlation function Ce{0) as a sum of a crossing-symmetric function C~^{6) 
and a crossing-antisymmetric function C^{6), 

CM = CM + c;(»), c|(») ^ c^W^CeI^-") . (4.8) 

Upon analytic continuation from the Euclidean to the Minkowskian theory, using 
Eq. (12.161) . one can show that they are related respectively to pomeron (i.e., C = +1) 
and odderon (i.e., C = —1) exchanges in the dipole-dipole scattering amplitude. A small 
but non-zero crossing-antisymmetric component is present in our data, thus signalling 
the presence of odderon contributions to the loop-loop correlation functions and in turn to 
the dipole-dipole scattering amplitudes. (Such contributions are absent in meson-meson 
elastic scattering [36l [23] , since in this case the relevant correlation function C|f^ is auto- 
matically crossing-symmetric, as already noticed at the end of Section 2). In Fig. [T7I we 
show the crossing-antisymmetric part of the ''zzz/zyy" data at d = 0, together with 
the corresponding SVM prediction. 



5. Conclusions and outlook 

In this paper the problem of high-energy hadron-hadron (dipole-dipole) scattering 
has been addressed (for the first time) from the point of view of lattice QCD. We have 
performed Monte Carlo numerical computations of Euclidean Wilson-loop correlation 
functions in SU (3) pure-gauge LGT. The energy dependence of soft scattering amplitudes 
at high energies is encoded in the dependence of these correlation functions on the relative 
angle 6 between the loops and can be reconstructed after the analytic continuation 6 — > 
—i log s/m^. 

In this paper we have focused on the study of the angular dependence of the Euclidean 
correlation function. An interesting and important feature that we have found is the pres- 
ence of an asymmetry with respect to ^ = 7r/2 in the plot of the Euclidean correlation 
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function against the relative angle: this (upon analytic continuation from Euclidean to 
Minkowskian theory) is a signal of the presence of C-odd contributions in dipole-dipole 
scattering, i.e., a signal of odderon exchange between the dipoles. Even though these C- 
odd contributions are averaged to zero in meson-meson scattering (at least in our model, 
as long as the squared meson wave functions satisfy some reasonable symmetry proper- 
ties in their dependence on the dipole orientations and on the longitudinal-momentum 
fractions), they might play a non-trivial role in more general hadron-hadron processes in 
which baryons and antibaryons are also involved. 

Although we cannot determine exactly the angular dependence from a finite set of 
numerical values, we can nevertheless compare the prediction of any given model with the 
lattice data by direct numerical comparison, if the model is quantitative, or by testing 
the given functional form with a best-fit to the data, if at least the 6'-dependence is 
known; in this way we can discriminate between various proposals and thus check the 
goodness of the approximations involved in the specific models. [One can of course also 
try to fit the data with some given arbitrary functions and then look at the results of the 
best-fits, but in accepting or rejecting a given function one must also take into account 
physical arguments, as two fitting functions can differ numerically by a small amount and 
nevertheless have different analytic structures that can result in completely different (and 
sometimes physically unacceptable!) high-energy behaviours after analytic continuation 
to Minkowski space-time.] 

The comparison of our data with the existing analytic calculations is not, generally 
speaking, fully satisfactory. 

• The SVM prediction (14. 3 p agrees with our lattice data in a few cases, at least in the 
shape and in the order of magnitude, but, in general, it is far from being satisfactory. 
More or less the same conclusion is reached if one instead performs a one-parameter 
best-fit with the given expression. This suggests that corrections to the SVM are 
in order, which could be relevant when deriving the high-energy behaviour of the 
scattering amplitudes, upon analytic continuation to Minkowski space-time. 

• We have then tried best-fits with a simple perturbative-like expression (14.41) and 
with the ILM expression (14.51) . The results are again not satisfactory. In particular, 
the ILM expression seems to be strongly disfavoured at c? = 0, while at (i = 2 it 
looks better than the SVM and perturbative-like expressions. (This suggests that 
while at d = the perturbative effects are dominant, ai d = 2 nonperturbative 
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effects are surely relevant and cannot be neglected.) By combining the two previous 
expressions into the ILMp expression fl4.6p . largely improved best-fits have been 
obtained. It would be interesting to also have a quantitative numerical prediction 
from the ILM and see how it compares with the lattice data. 

• Finally, we have tried a best-fit with the AdS / CFT expression fl4.7p , although there 
is no reason for this formula to apply in our case, since it was originally derived 
for the A/" = 4 SYM theory at strong coupling and large impact parameter. Taking 
into account that this is a three-parameter best-fit, even this one is not satisfac- 
tory, especially at d = 0. Best-fits with QCD-inspired expressions with only two 
parameters, like, e.g., the ILMp expression (14.61) [or some appropriate modification 
of the SVM expression fl4.3l) ] give smaller Xd.o.f.- 

Although the AdS/CFT expression (14.71) . as said above, is not expected to describe real 
QCD, it nevertheless shows how a non-trivial high-energy behaviour could emerge from 
a simple analytic dependence on the angle 6. However, in this case, after the analytic 
continuation into Minkowski space-time, it is not possible to pass to the high-energy 
limit under the integral sign, as the integrand is an oscillating function of the energy, 
and one should carry over the remaining integrals first. As we have already pointed 
out, the integration over the distance between the loops should be the relevant one and, 
depending on the detailed form of the various coefficient functions, a variety of behaviours 
could emerge. It seems then worth investigating further the dependence of the correlation 
functions on the relative distance between the loops, as well as on the dependence on 
the relative angle, as they could combine in a non-trivial way to determine the high- 
energy behaviour of meson-meson total cross sections. Moreover, as already recalled in 
the previous Section [see Eq. (14. 2p ]. the study of the transverse-distance dependence of 
the Euclidean correlation function Qe 6 = 0,7t would allow to determine the static 
dipole-dipole potential. These and other related issues will be addressed in future works. 

In conclusion, the fact that the existing models are not able to fully explain the lattice 
data, which inside the errors represent the true QCD expectation, is a motivation to 
further investigate the problem of soft high-energy scattering, both on the numerical and 
on the analytical side. We hope that the interplay of numerical simulations in LGT and 
analytic nonperturbative calculations can lead to a deeper understanding of the long- 
standing problem of high-energy scattering in strong interactions. 
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FIGURE CAPTIONS 



Fig. 1 The space-time configuration of tlie two Wilson loops Wi and W2 entering in the ex- 
pression for the dipole-dipole elastic scattering amplitude in the high-energy limit. 

Fig. 2 The analyticity domains of the function in the complex variable 9 and of the 
function Qm in the complex variable x- 

Fig. 3 The minimal-distance prescription for a line with slope tan^ = 1/2. 

Fig. 4 The relevant Wilson-loop configuration. Using the 0(4) invariance of the Euclidean 
theory we have put piE parallel to the xei axis. 

Fig. 5 Loop configuration in the transverse plane. In the "af e" case the link orientation is 
not shown as it is averaged over. 

Fig. 6 Dependence of Qe on the length Li = L2 = L (in lattice units) of the loops at 
9 = 90° for rf = 0,1,2. 

Fig. 7 Angular dependence oi Qe for various lengths of the loops, for = 0. The numbers 
in the key refer to the on-axis loop length and to the off-axis loop set respectively, 
as explained in the text and in Table 1. (Different data sets have been slightly 
shifted horizontally for clarity.) 

Fig. 8 Angular dependence oi Qe for various lengths of the loops, for d = 1. The notation 
is the same as in Fig. [71 

Fig. 9 Angular dependence oi Qe for various lengths of the loops, for d = 2. The notation 
is the same as in Fig. [3 

Fig. 10 Dependence of CE^d) (on a logarithmic scale) on the distance d (in lattice units) at 
9 = 45° and 9 = 90°. 

Fig. 11 Comparison of the lattice data to the SVM prediction (14.31) with i^svM calculated 
according to Ref. [TS] (solid line) and to the one-parameter (-R'svm) best -fit (for the 
"2:2:2" and "21/1/' cases only) with the SVM expression (14.31) (dotted line) at c? = 0. 

Fig. 12 The same comparison as in Fig. [11], but at d = 1. 
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Fig. 13 The same comparison as in Fig. [TT], but at d = 2. 

Fig. 14 Comparison of lattice data to best-fits with the perturbative-like expression fl4.4p 
(solid line), the ILM expression (14.51) (dotted line) and the ILMp expression (14. 6 p 
(dashed line) at d = 0. 

Fig. 15 The same comparison as in Fig. [HI but at = 1. 

Fig. 16 Comparison of lattice data to a best-fit with the AdS/CFT expression (14. 7p for 
various cases. 

Fig. 17 The crossing-antisymmetric component Cg, as defined in Eq. (14. 8p . for the "zzz/ zyy" 
case at d = and the corresponding prediction using the SVM expression (14.31) (solid 
line) . 
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Figure 1: The space-time configuration of tlie two Wilson loops Wi and VV2 entering in 
the expression for the dipole-dipole elastic scattering amplitude in the high-energy limit. 




Figure 2: The analyticity domains of the function in the complex variable 6 and of 
the function Q m in the complex variable x- 



Figure 3: The minimal-distance prescription for a line with slope tan^ 
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Figure 4: The relevant Wilson-loop configuration. Using the 0(4) invariance of the 
Euclidean theory we have put piE parallel to the xei axis. 
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Figure 5: Loop configuration in the transverse plane. In the "ave" case the hnk orientation 
is not shown as it is averaged over. 
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Figure 6: Dependence of Qe on the length Li = L2 = L (in lattice units) of the loops at 
e = 90° for rf = 0,1,2. 
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Figure 7: Angular dependence oi Qe for various lengths of the loops, for d = 0. The 
numbers in the key refer to the on-axis loop length and to the off-axis loop set respectively, 
as explained in the text and in Table [H (Different data sets have been slightly shifted 
horizontally for clarity.) 
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Figure 8: Angular dependence of Qe for various lengths of the loops, for (i = 1. The 
notation is the same as in Fig. [71 



0.005 



0.004 



0.003 



0.002 



0.001 



, 1 

n 


1 


1 1 


1 1 


1 










, g 














□ 






J I* . 




1 




: f 












1 


1 


1 1 


1 

1 1 


1 



0.005 



0.004 



g 0.003 
0.002 



0.001 



0.005 



0.004 



0.003 



0.002 



0.001 



6,1 I— H- 
6 2 ---X- 
8 1 

8 2 - ■□■ 



30 


50 




70 


90 

en 
zyy 


110 


130 


150 


1 


1 




1 


1 


1 


1 


1 

□ 


-, n 
















X ■ 














X 


- it 


: o 

Xi' 










, n 
1 


it _ 




J* 


: □ 

X:, 

St 


1 




1 


1 




30 


50 




70 


90 

en 

ave 


110 


130 


150 



6,1 I— H- 
6 2 ■ X 
8 1 '--51^- 
8 2 - d- 




Figure 9: Angular dependence of Qe for various lengths of the loops, for d = 2. The 
notation is the same as in Fig. [71 
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Figure 10: Dependence of C£;((i) (on a logarithmic scale) on the distance d (in lattice 
units) at ^ = 45° and ^ = 90° . 
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Figure 11: Comparison of the lattice data to the SVM prediction (14.31) with K^ym cal- 
culated according to Ref. [18] (solid line) and to the one-parameter (-R'svm) best -fit (for 
the and "2?/?/' cases only) with the SVM expression (14.31) (dotted line) at c? = 0. 





Figure 13: The same comparison as in Fig. [TT|, but at ci = 2. 




Figure 14: Comparison of lattice data to best-fits with the perturbative-hke expression 
(14.41) (sohd hne), the ILM expression (14.51) (dotted hne) and the ILMp expression (14. 6 p 
(dashed hne) at c? = 0. 




Figure 15: The same comparison as in Fig. [HI but at ci = 1. 
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Figure 16: Comparison of lattice data to a best-fit with the AdS/CFT expression (14. 7p 
for various cases. 
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Figure 17: The crossing-antisymmetric component C^, as defined in Eq. (14. 8p . for the 
"zzz/ zytf^ case ai d = and the corresponding prediction using the SVM expression (14. Sp 
(solid line). 



